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Rickettsia aeschlimannii is a tick-associated human pathogen. We report here the draft genome of R. aeschlimannii strain 
MC16, isolated from Hyalomma marginatum marginatum ticks collected in Morocco. 
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Rickettsia spp. are strictly intracellular bacteria that are associ- 
fl ated with arthropods. Some of them cause human infections 
that are currently considered emerging and have been identified 
worldwide (1). Rickettsia aeschlimannii is a member of spotted 
fever group rickettsiae that was first isolated from Hyalomma mar- 
ginatum marginatum ticks collected from camels in Morocco in 
1992 (2). The first human infection was reported in 2002 in a 
French patient who had traveled to Morocco (3). R. aeschlimannii 
has since been reported to cause infections in South Africa (4), 
Tunisia (5), Algeria (6), and Greece (7). Here, we briefly describe 
the draft genome of R. aeschlimannii strain MC16 T . 

The MiSeq sequencer (Illumina, San Diego, CA) was used to 
sequence the genomic DNA of R. aeschlimannii MC16 (depos- 
ited in the Collection de Souches de FUnite des Rickettsies 
under reference CSUR R8) using a mate-pair library. Quality 
trimming and de novo assembly of the reads were performed 
using the CLC Genomics Workbench version 6.0.1 (CLC bio, 
Aarhus, Denmark). The resulting contigs were reordered in 
Mauve version 2.3.1 (8) using Rickettsia montanensis strain 
OSU 85-930 as a reference genome (accession no. CP003340.1). Po- 
tential coding sequences (CDSs) were predicted using AMIGene 
(9), and the assignment of protein functions was performed by 
searching against the RickBase (10), GenBank, and Pfam (11) da- 
tabases using BLASTp (12), while rRNAs, tRNAs, and other RNAs 
were identified using BLASTn, tRNAscan-SE version 1.21 (13), 
and RNAmmer 1.2 (14), respectively. Orthologous genes between 
R. aeschlimannii and R. montanensis were identified using Or- 
thoMCL (15), with a BLASTp E value cutoff of 1 X 10" 5 and the 
default Markov cluster algorithm (MCL) inflation parameter of 
1.5. 

The draft genome of R. aeschlimannii MC16 consists of 16 
contigs ranging in size from 314 to 238,303 bases, resulting in a 
total genome of 1,312,196 nucleotides, with an average genome 
coverage of 145-fold and a G+C content of 32.2%. Two plasmids 
of 14,948 bp and 26,886 bp, respectively, were identified and ex- 
hibited identity matches of 96% (97% coverage; E value, 0.0) 
to Rickettsia rhipicephali strain 3-7-female6-CWPP plasmid 
pMCC_l (accession no. CP003343.1) and 96% (27% coverage; 
E value, 0.0) to Rickettsia monacensis strain IrR/Munich plasmid 



pRM (accession no. EF564599.1), respectively. The chromosome 
contains 1,790 CDSs and, like other rickettsial genomes, 3 non- 
contiguous rRNAs (5S, 16S, and 23S rRNAs), 33 tRNAs, and 3 
other RNAs. In addition, the two plasmids contain 23 and 44 
CDSs but no RNAs. 

The R. aeschlimannii chromosome exhibits a high level of syn- 
teny to R. montanensis but differs from it by three inversions of 
66,949 bp, 3,404 bp, and 53,495 bp and the absence of the genes 
encoding the integration host factor beta-subunit (himD), large 
extracellular a-helical protein, 3-hydroxyacyl-coenzyme A (CoA) 
dehydrogenase {fadB), and competence protein F {comFl and 
comF2). 

Nucleotide sequence accession numbers. This whole-genome 
shotgun project has been deposited at DDBJ/EMBL/GenBank un- 
der accession no. CCEA01000001 to CCEA01000014 (BioProject 
PRJEB6087), and the two plasmid sequences have been deposited 
under accession no. CCEB01000001 and CCEB01000002. 

ACKNOWLEDGMENT 

This research was funded by the Mediterranee-lnfection foundation. 
REFERENCES 

1. Parola P, Paddock CD, Socolovschi C, Labruna MB, Mediannikov O, 
Kernif T, Abdad MY, Stenos J, Bitam I, Fournier PE, Raoult D. 2013. 
Update on tick-borne rickettsioses around the world: a geographic ap- 
proach. Clin. Microbiol. Rev. 26:657-702. http://dx.doi.org/10.1128/ 
CMR.00032-13. 

2. Beati L, Meskini M, Thiers B, Raoult D. 1997. Rickettsia aeschlimannii 
sp. nov., a new spotted fever group rickettsia associated with Hyalomma 
marginatum ticks. Int. J. Syst. Bacteriol. 47:548-554. http://dx.doi.org/ 
10. 1099/00207713-47-2-548. 

3. Raoult D, Fournier PE, Abboud P, Caron F. 2002. First documented 
human Rickettsia aeschlimannii infection. Emerg. Infect. Dis. 8:748-749. 
http://dx.doi.org/10.3201/eid0807.010480. 

4. Pretorius AM, Birtles RJ. 2002. Rickettsia aeschlimannii: a new patho- 
genic spotted fever group rickettsia, South Africa. Emerg. Infect. Dis. 
8:874. http://dx.doi.org/10.3201/eid0808.020207. 

5. Znazen A, Rolain JM, Hammami A, Jemaa MB, Raoult D. 2006. 
Rickettsia felts infection, Tunisia. Emerg. Infect. Dis. 12:138-140. http:// 
dx.doi.org/10.3201/eidl201.050876. 

6. Mokrani N, Parola P, Tebbal S, Dalichaouche M, Aouati A, Raoult D. 
2008. Rickettsia aeschlimannii infection, Algeria. Emerg. Infect. Dis. 14: 
1814-1815. http://dx.doi.org/10.3201/eidl411.071221. 



July/August 2014 Volume 2 Issue 4 e00666-14 



Genome Announcements 



genomea.asm.org 1 



Sentausa et al. 



7. Germanakis A, Chochlakis D, Angelakis E, Tselentis Y, Psaroulaki A. 
2013. Rickettsia aeschlimannii infection in a man, Greece. Emerg. Infect. 
Dis. 19:1176-1177. http://dx.doi.org/10.3201/eidl907.130232. 

8. Rissman AI, Mau B, Biehl BS, Darling AE, Glasner JD, Perna NT. 2009. 
Reordering contigs of draft genomes using the Mauve aligner. Bioinfor- 
matics 25:2071-2073. http://dx.doi.org/10.1093/bioinformatics/btp356. 

9. Bocs S, Cruveffler S, Vallenet D, Nuel G, Medigue C. 2003. AMIGene: 
annotation of Microbial genes. Nucleic Acids Res. 31:3723-3726. http:// 
dx.doi.org/10.1093/nar/gkg590. 

10. Blanc G, Ogata H, Robert C, Audic S, Suhre K, Vestris G, Claverie JM, 
Raoult D. 2007. Reductive genome evolution from the mother of Rickett- 
sia. PLoS Genet. 3:el4. http://dx.doi.org/10.1371/journal.pgen.0030014. 

11. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang 
N, Forslund K, Ceric G, Clements J, Heger A, Holm L, Sonnhammer 
ELL, Eddy SR, Bateman A, Finn RD. 2012. The Pfam protein families 



database. Nucleic Acids Res. 40:D290-D301. http://dx.doi.org/10.1093/ 
nar/gkrl065. 

12. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, 
Lipman DJ. 1997. Gapped BLAST and PSI-BLAST: a new generation of 
protein database search programs. Nucleic Acids Res. 25:3389-3402. 
http://dx.doi.org/10.1093/nar/25.17.3389. 

13. Lowe TM, Eddy SR. 1997. tRNAscan-SE: a program for improved detec- 
tion of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25: 
955-964. http://dx.doi.Org/10.1093/nar/25.5.0955. 

14. Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. 
2007. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. 
Nucleic Acids Res. 35:3100-3108. http://dx.doi.org/10.1093/nar/gkml60. 

15. Li L, Stoeckert CJ, Roos DS. 2003. OrthoMCL: identification of ortholog 
groups for eukaryotic genomes. Genome Res. 13:2178—2189. http:// 
dx.doi.org/10.1101/gr.1224503. 



2 genomea.asm.org 



Genome Announcements 



July/August 2014 Volume 2 Issue 4 e00666-14 



